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ABSTRACT 

We report two-dimensional spectroastrometry of Br7 emission of TW Hya to study the kinematics 
of the ionized gas in the star-disk interface region. The spectroastrometry with the integral field 
spectrograph SINFONI at the Very Large Telescope is sensitive to the positional offset of the line 
emission down to the physical scale of the stellar diameter (~0.01 AU). The centroid of Br7 emission 
is displaced to the north with respect to the central star at the blue side of the emission line, and 
to the south at the red side. The major axis of the centroid motion is P.A.= —20°, which is nearly 
equal to the major axis of the protoplanetary disk projected on the sky, previously reported by CO 
sub millimeter spectroscopy (P.A.= —27°) The line-of-sight motion of the Br7 emission, in which the 
northern side of the disk is approaching toward us, is also consistent with the direction of the disk 
rotation known from the CO observation. The agreement implies that the kinematics of Br7 emission 
is accounted for by the ionized gas in the inner edge of the disk. A simple modeling of the astrometry, 
however, indicates that the accretion inflow similarly well reproduces the centroid displacements of 
Br7, but only if the position angles of the centroid motion and the projected disk ellipse is a chance 
coincidence. No clear evidence of disk wind is found. 

Subject headings: circumstellar matter — planetary systems: protoplanetary disks — stars: activity 
— stars: formation — stars: individual (TW Hya) — stars: pre-main sequence — 
techniques: high angular resolution 



1. INTRODUCTION 

It is being commonly accepted that a protoplanetary 
disk dissipates inside out. Although the direct trigger 
is still not clear, a large number of spectral energy dis- 
tributions collected by the Spitzer Space Telescope (e.g., 
Currie et al. 2008; Sicilia-Aguilar et al. 2008), as well as 
the direct imaging of the central cavity by the submil- 
limeter/millimeter interfcrometric arrays (submillimeter 
array (SMA): Brown et al. 2008; Plateau de Bure In- 
terferometer: Pietu et al. 2006), give clear observational 
supports that many disks dissipate at the central part 
first. The innermost region of a disk (< 1 AU) is, there- 
fore, the critical site to study how the disk dissipation 
starts. The structure of a disk at that radius is, however, 
expected to be complicated (see Dullemond & Monnier 
2010, for a recent review). Dust grains evaporate in 
close proximity to the star. Hot gas may remain inside 
the inner truncation of the dust disk, but in ongoing dis- 
sipation either by accretion onto the star, disk wind/jet, 
or by a forming planet that sweeps out the disk. 

TW Hya belongs to a small group of stars that the disk 
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dissipation is currently in progress (Alexander 2007, for 
a review of transition objects). Close vicinity of TW Hya 
to the solar system (50-60 pc; Kastner et al. 1997; Hoff, 
Henning, & Pfau 1998; Mamajek 2005) makes it an ex- 
ceptional laboratory to look into the central few AU 
of a disk (we will use d =51 pc from Mamajek (2005) 
throughout this paper). Missing inner disk of TW Hya 
was first recognized by Calvet et al. (2002) through the 
lack of near-infrared excess, which attests the absence 
of hot dust close to the central star. A comprehensive 
multiwavelength study of the mass distribution of the 
disk by Gorti et al. (2011) confirmed that the inner disk 
(< 3.5 AU) is indeed depleted by one to two orders of 
magnitude than that is extrapolated from the surface 
density of the disk immediately outside. Interferometric 
observations also lend support to the clearing of the inner 
disk, although the exact location of the inner truncation 
is under debate (0.4-4 AU; Eisner et al. 2006; Ratzka et 
al. 2007; Hughes et al. 2007). 

Hydrogen recombination lines are widely used to study 
the hot plasma near the central star, in particular in 
the accretion columns channeled by the magnetic fields 
(Uchida & Shibata 1984; Camenzind 1990; Konigl 1991). 
However, the quantitative modeling of the line emission 
is still a challenge (Romanova et al. 2004; Gregory et al. 
2006; Kurosawa, Harries, & Symington 2006; Kurosawa 
et al. 2008), and the origin of the emission line is not 
completely understood, either (e.g., Kraus et al. 2008). 
The goal of this paper is to study the dynamics of the ion- 
ized gas inside the inner hole of a transition disk with a 
high angular precision afforded by two-dimensional (2D) 
spectroastrometry (e.g., Garcia et al. 1999; Davies et 
al. 2010). Spectroastrometry with long slit spectrograph 
has been used for some time to study the kinematics of 
the emission lines of pre-main-sequence stars on the or- 
der of 1 AU (e.g., Solf & Bohm 1993; Bailey 1998). We 
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Fig. 1. — Top: the centroids of the PSF at the wavelengths near 3vy emission. North is up, and east is to the left. The centroids recorded 
with the detector column aligned to P.A.=0° are shown in the left panel, and to P.A.=90° in the right panel. The PSF centroids in the 
Br7 emission are shown with squares, and those of the continuum wavelengths are shown with crosses. Both symbols are color-coded from 
blue to red with the wavelength. The drifts of the PSF at the continuum wavelengths are fitted by linear functions of the wavelength, and 
are shown in dotted lines. Bottom: line profiles of the Br7 emission. The crosses and the squares are from the observations, color-coded 
in the same way with the panels above. The red lines are Lorentzian functions that best fit the observed data points. 



will employ an integral field spectrograph assisted by an 
adaptive optics system to achieve the astrometry of the 
physical scale of a stellar diameter in order to deal with 
the ionized gas in the star-disk interface region. 

2. OBSERVATION 

The observations were performed by the integral field 
spectrograph SINFONI (Eisenhauer et al. 2003) at the 
Very Large Telescope (VLT) with the adaptive optics sys- 
tem MACAO (Bonnet et al. 2004). An image on the focal 
plane (field of view 0'.'8 x 0'.'8) was sliced into 32 reflective 
slitlets, each providing a medium resolution spectrum in 
the K band (R =4000). The pixel scale is 12.5 mas on the 
detector, although the spatial sampling perpendicular to 
the length of the slitlets is limited by the resolution of the 
image sheer (25 mas). The observations were carried out 
with two instrumental angles, one with the detector col- 
umn pointing to the north (P.A.=0°; performed on UT 
2006 January 19), and the other to the east (P.A.=90°; 
on UT 2006 January 28). 

3. DATA REDUCTION 

3.1. Spectroscopy 

The SINFONI data-reduction pipeline provides data 
cubes that the monochromatic images are stacked along 
the wavelength. One-dimensional (ID) spectra were ex- 
tracted from the data cubes by integrating the flux in- 
side the aperture of 0'.'2. The spectroscopic standard 



star HD 37702 (B8 V) was used to remove the telluric 
lines. The photospheric absorption line of HD 37702 at 
Br7 was removed before dividing the object spectra by 
subtracting a Lorentzian function fitting the absorption 
line profile. The wavelength calibration was performed 
by cross-correlating the object spectra with the atmo- 
spheric transmission curve calculated by ATRAN (Lord 
1992). The error of the wavelength calibration was eval- 
uated against the photospheric absorption lines at 2.10- 
2.21 yum (Al I, Mg I, Na I), and found to be accurate to 
20 km s . The observed line velocity was converted to 
the velocity with respect to the local standard of rest, 
using rv package on IRAF. 8 The equivalent width of the 
line emission was reduced to ~ 30% of that on 2006 Jan- 
uary 19 (P.A.=0°), when we observed the object again 
on 2006 January 28 with the instrumental position angle 
90°. The variability of Br7 emission of a similar am- 
plitude has been reported by Eisner et al. (2011) for 
TW Hya. 

3.2. Drift Correction 

The point-spread function (PSF) of the star near Br7 
emission was fit with a 2D Gaussian function at each 
wavelength to measure the PSF centroid. The positions 
of the PSF centroid have shown an overall drift along 

8 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 2. — Same as in Figure 1, but after the correction for the PSF drift. The offsets of the PSF centroids are converted to the physical 
scale, assuming the distance to TW Hya is 51 pc (Mamajck 2005). The major axis of the locus is oriented to P.A.=— 20° (left). 



P. A. =71° with the wavelength in the case of the instru- 
mental angle aligned to 0° (crosses in Figure 1, left). This 
is apparently an instrumental effect, left uncorrected by 
the image alignment process in the pipeline, as the drift 
orientation changes by about the same amount when the 
instrument is rotated by 90° (the drift P.A.= -10°, Fig- 
ure 1, right). The overall motion of the PSF drift was 
fit by a linear function of the wavelength and subtracted 
from the measurements. The centroid position at Br7 
emission after the correction is shown in Figure 2. The 
motion of the emission centroid of B1-7 was unfortunately 
aligned to the orientation of the PSF drift in the data 
with the instrumental angle set to P.A.=90°. Although 
the results of the two observations performed with differ- 
ent rotation angles are qualitatively consistent, we will 
use the data from P.A.=0° mostly in the following dis- 
cussion, as it is less affected by the correction of the PSF 
drift. The locus of the centroid motion is an elongated 
ellipse with the major axis aligned to P.A.= —20° in the 
data obtained with P.A.=0°. The total extent of the loop 
is about 0.05 pixel, or 0.04 AU (8i?*) from the central 
star at the distance of TW Hya (i?*=0.96 Rq; Yang, 
Johns-Krull, & Valenti 2005). 

3.3. Astrometric Accuracy 

The smallest centroid offset that we can detect in 
the observation was estimated from the dispersion of 
the measured centroid positions at the continuum wave- 
lengths. The standard deviation of the amplitudes of the 
centroid positions at the continuum (Figure 2) is 5 x 10 -3 
pixel, which is 0.003 AU at the distance of TW Hya, 
or less than one stellar radius. We will check below if 
the high astrometric accuracy achieved is consistent with 
the simple photon statistics. In the case that the error 
distribution is canonical, the accuracy of the astrome- 
try improves from the size of PSF proportionally to the 
squared root of the total photons detected (e.g., Bai- 
ley 1998). The size of the PSF during the observation 
of TW Hya measured in the median frame of the im- 
age stack was 65 mas in the full width at half-maximum 
(FWHM) . The integrated flux of the star in a raw frame 
was ~lx 10 5 ADU for 30 s of exposure time at the con- 



tinuum wavelengths near Br7 emission. The total expo- 
sure time was 360 s for the object for one instrumental 
rotation. The photon count of the fully processed im- 
age we used for the astrometry is therefore ~ 2.9 x 10 6 
per spectral unit, after converting ADU to photons by 
the detector gain 2.4 ADU/electron. The nominal sta- 
tistical error of the astrometry is therefore 0.038 mas or 
3.1 xl0~ 3 pixel. The astrometric accuracy achieved is 
about 1 x 10~ 2 pixel in terms of the FWHM, therefore 
consistent with the photon statistics, though it is a few 
times worse than the ideal case. 

We can calculate a formal error of the center of gravity 
for more general cases, using the realistic instrument pa- 
rameters of SINFONI. The flux-weighted centroid of an 
emission source is given by 



z(X) = 



i 



where Xi is the ID coordinate of the detector either row 
or column, and yi is the photon count of the pixel Xi. 
The formal error of the centroid position a z is given by 
the pixel counts and their errors a yi , 



E 



i i 

Q» 2 



We calculated a Vi as a squared sum of the readout noise 
and the shot noise of the pixel, 

a li = ( a r + Vi) ■ n h, 

where rify is the number of frames averaged. The read- 
out noise er r is about 10 electron for a single readout 
for SINFONI. The pixel count yi is evaluated with ac- 
tual observed data in the unit of electrons, and before 
the subtraction of the sky emission; although the sky 
background is negligible at the wavelength concerned (0- 
10 ADU). The formal error er 2 is a function of the size 
of the area in which the summation is calculated and 
is smallest when the minimum number of the pixels is 
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Fig. 3. — Comparison of the astrometry measured by Gaussian fitting (left) to that by the center of gravity (right). The center of gravity 
(flux- weighted center of the emission) is calculated within 8x8 pixel windows centered at the peak of the PSF. Two measurements are 
consistent, attesting the positive astrometric signal at Br7 emission is not an artifact of the measuring scheme. 



involved. a z is on the order of 0.01 pixel from 4x4 to 
16x16 pixel windows, which is in reasonable agreement 
with 0.005 pixel that is actually achieved. The center of 
gravity measured in 8 x 8 pixel window is compared with 
the measurements of the centroid by 2D Gaussian fitting 
and found to be consistent (Figure 3). 

3.4. Test Against Artifact 

A detection of a spectroastrometric signal is often 
tested with photospheric emission lines, in which little 
or no offset of the emission centroid is expected. Un- 
fortunately, we cannot afford such a test, because Br7 
emission is the only significant emission line in our spec- 
tral coverage. Instead, we measured the astrometry of 
PSF standard stars at the wavelength of Br7 line with 
the same method that we used for TW Hya in order 
to show that the offset of the centroid detected is not 
an artifact of the measuring scheme. The PSF stan- 
dard stars StKM 1-440 and CE 284 were observed on UT 
2006 January 24 and 28, respectively. They are fainter 
than TW Hya (#=7.297 mag; StKM 1-440 #=8.98 mag; 
CE 284 #=9.39 mag) and were observed with less tele- 
scope time. The astrometric accuracy achieved was 
therefore a few times worse, as the photon statistics are 
not as good as TW Hya. Nevertheless, the astrometric 
accuracy is good enough to negate the possible offset of 
the centroid at the wavelength of Br7, in the sense that 
it should have been detected, if the astrometric signal in 
TW Hya were an artifact and were present in the data 
of StKM 1-440 and CE 284 with the same amplitudes 
(Figure 4). 

The advantage of 2D spectroastrometry is not only a 
better spatial coverage and better photon statistics, but 
also the robustness against the artifacts. The most of the 
artifacts that spectroastrometry suffers from, come from 
the use of a slit and vignetting of the source, combined 
with uneven illumination of a slit and a distorted PSF 
(Brannigan et al. 2006). In 2D spectroastrometry, we do 
not use a slit but slitlets, where the problem of vignetting 
is expected to be much smaller. 

3.5. Convolution with Continuum Emission 



Note that the displacement of the centroid of Br7 emis- 
sion we measured above is a lower limit of the real offset 
of the emission line, because the PSF of the line emission 
is convolved with that of the continuum emission. The 
convolution shifts the line emission centroid toward the 
continuum emission source, if it is measured by fitting 
a Gaussian function. If we assume that the amount of 
the shift is proportional to the flux ratio, the convolution 
moves the observed centroid to 

r ° bS = F ,+Fr ' 

± cont ~ ± line 

where ri; no is the intrinsic location of the centroid of 
the Br7 emission before the convolution, and Fcont and 
inline a re the flux of the continuum and the continuum- 
subtracted line emission, respectively. Effect of such con- 
volution can be seen in Figure 1, where the amplitude of 
the centroid motion is smaller in the data with the instru- 
mental angle set to P.A =90° than that to P.A.=0°. This 
is reasonable, since the line intensity with respect to the 
continuum emission is three times smaller in P.A. =90°. 
As we are only interested in the location of the line cen- 
troid relative to that of the continuum emission which 
is supposed to be the position of the star, we can set 
r CO nt = and r linc is given by 

^"linc — ^obs(l ~\~ -^cont /-^line) • 

Such a simple correction, though frequently used, is prob- 
lematic, because the correction factor 1 + F cont /Fy mc is 
large at the line wings where the line emission is weak 
compared to the continuum emission (-Fiine/-Fcont ~ 0); 
and the astrometric signal, as well as its error, diverge 
into infinity at the continuum wavelengths. 9 To avoid 
that, one could apply the correction only at the wave- 
lengths where the line emission is detected above a cer- 
tain significance (-Fiin e /-Fcont 3> 0). However, such a se- 
lective correction has a pitfall to distort the geometry 
of astrometry by exaggerating astrometric signals at the 
border of the wavelength interval that the corrections 
were applied. 

9 In our case, the line-to-continuum ratio at the line center is 0.75 
(P.A.=0°) and 0.25 (P.A. =90°), therefore the correction factor is 
at minimum 1 + F con t /Fn nc > 2.3 and > 5. The correction factor 
at the line wings is much larger than the line center. 
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Fig. 4.— Comparison of the astrometry of TW Hya with P.A.=0° to that of two PSF standard stars StKM 1-440 and CE 284 at the 
same wavelength with Br7 emission. Note that the astrometric accuracy of the PSF stars is not as good as TW Hya, because the photon 
statistic is worse. Nevertheless, it is reasonably clear that no positive astrometric signal is detected in either of the PSF standard stars, if 
the amplitude of the possible false signal is as large as TW Hya. 



TABLE 1 

Geometrical Parameters of TW Hya Stellar-Disk System 
Assumed Throughout the paper. 



Parameter Values Reference 

Stellar mass 0.8 M Q Calvet et al. (2002) 

Stellar radius 0.96 Rq Yang, Johns-Krull, & Valenti (2005) 

Co-rotational radius 6.3 R t Johns-Krull & Valenti (2001) 

Disk P.A. a —20° Present observation 

Disk inclination angle 7° Qi et al. (2004, 2006) 

Distance 51 pc Mamajek (2005) 



is almost face-on (i <4°), according to the imaging in 
the scattered light by Hubble Space Teles cope/ WFPC2 
and NICMOS (Krist et al. 2000; Weinberger et al. 2002). 
Pontoppidan et al. (2008) performed ID spectroastrome- 
try on CO fundamental line and proposed similarly small 
inclination angle (i = 4°±1°). We use i = 7° reported in 
the submillimeter CO observation by the SMA (Qi et al. 
2004, 2006) in the following models, which is constrained 
in spatially and spectroscopically resolved manner. The 
position angle of the disk major axis is set to P.A.=— 20°, 
equal to the position angle of the loop observed in the 
present study. The stellar radius (R* — 0.96 Rq) and the 
mass of the star (M* = 0.8 Mq) are taken from Yang, 
Johns-Krull, & Valenti (2005) and Calvet et al. (2002), 
respectively. The disk is assumed to be truncated inward 
at the co-rotational radius R c = 6.3 i?* (Johns-Krull & 
Valenti 2001). The stellar and disk parameters used in 
the models are summarized in Table 1. 

The centroid of Br7 emission in this model is displaced 
to the north at the blueshifted wing, and moves across 
the disk with the wavelength along the position angle 
perpendicular to the disk rotation axis. This is qual- 
itatively consistent with the observation (Figure 6(a)). 
However, the apparent radial velocity of the gas is too 
slow to match the observation. Keplerian velocity at the 
co-rotation radius is v sini < 19 km s _1 . After convolved 
with the thermal broadening (v =9 km s _1 for 10 4 K) 
and the instrumental resolution (v —75 km s -1 ), the line 
width of Br7 emission is still too narrow by itself to ex- 
plain the observed line width of 150 km s _1 in FWHM 
(Figure 6(a)). 

4.2. Accretion Inflow 

In the second model, the ionized gas is in a pair of 
accretion columns starting on a sphere at the distance of 
the co-rotation radius. The gas free falls onto the star 
in straight lines symmetrically placed about the central 



a The position angle of the disk plane projected on the sky. 

For these reasons, we choose not to apply the correc- 
tion for the continuum convolution with the understand- 
ing that we will refrain from the quantitative discussion 
concerning the physical scale and the complex aspects of 
the geometry of the ionized gas. We instead focus on the 
position angle of the displacement, which is more critical 
for the interpretation of the astrometric signal and is less 
affected by the convolution with the continuum emission. 

4. DISCUSSION 
4.1. Rotating Disk 

We use three simple geometric models of the ionized 
gas to compare to the observed centroid motion of B1-7 
emission: the gas at the inner rim of a rotating disk, in 
an accretion flow, and in a disk wind. 

In the rotating disk model, the plasma is at the inner 
edge of a thin disk (Figure 5(a)). The gas is assumed to 
be in pure Keplerian rotation with the apparent radial 
velocity given by v e = GM„ / R c cos 6 sin (— i) , where 
i and 9 are the inclination and the azimuthal angle on 
the disk plane, respectively. The geometry of the disk 
around TW Hya is reasonably well known. The disk 
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Fig. 5. — Schematics of the geometry assumed in the calculations of the model astrometry shown in Figure 6. The figures on the bottoms 
are the model disks on the top row, but in projection on the sky with the proper position angle of the disk plane. 



star (Figure 5(b)). The radial velocity component of the 

free- falling gas is v r = y/2GM*(l/r — 1/i?*) cos (i + b) 
at the distance of r with the impact angle b from the disk 
plane. The misalignment between the accretion columns 
and the disk rotation axis is assumed to be small, with 
b = 70°, so that the maximum radial velocity matches 
to the broad wings of Br7 emission (±200 km s _1 ). The 
azimuthal angle of the inflow projected on the disk is 
set to —20°, so that the apparent position angle of the 
astrometry is consistent with the observations. 

Br7 emission appears most blueshifted at the central 
star, as the gas in free-fall is accelerated until it hits the 
stellar surface. The centroid of the line emission moves 
outward with the wavelength along the accretion column 
up to the co-rotational radius. The line emission then 
jumps to the other end of the co-rotational circle to the 
near-side accretion column at the red half of the emis- 
sion line (Figure 6(b)). The apparent motion of the line 
centroid is qualitatively consistent with the observation. 
However, it remains unexplained how the position angle 
of the accretion inflow coincides with that of the disk 
major axis at the time of the observation. The disk of 
TW Hya is almost pole on. If the accretion columns are 
a single pair of stream lines tied to the stellar surface 
by the magnetic field, the position angle projected on 
the sky can be of any direction at a given moment, as it 
rotates with the stellar rotation. 

4.3. Disk Wind 

In the third model, the disk wind launches outward 
at the co-rotation radius with the launching angle a = 
30° (Figure 5(c)). The launching velocity is arbitrar- 
ily set to vq =100 km s _1 . The radial component of 
the disk wind at the disk azimuthal angle 9 is vg = 
— sin a cos i — vq cos a sin (i) sin (0) at the near-side of 
the disk. 

The centroid of B1-7 emission at the bluest wing ap- 
pears at the intersection of the co-rotational circle and 
the disk rotation axis projected on the sky (Figure 6(c)). 
The centroids of the line emission moves along the disk 
rotational axis with the wavelength from one end to the 
other end of the co-rotation circle. The orientation of 



the motion is perpendicular to what is observed. Disk 
winds should only have a limited contribution to B1-7 
kinematics on the physical scale we discuss here. 

5. CONCLUDING REMARK 

We conclude that the observed astrometry of the Br7 
emission can be accounted for either by disk rotation 
or accretion inflow. The disk rotation, however, does 
not explain the line width of Br7 emission by itself 
(150 km s _1 in FWHM), and the accretion inflow re- 
quires arbitrary assumption that the position angle of 
the accretion columns are aligned to the major axis of 
the disk ellipse at the time of the observation. No clear 
evidence for disk winds is found. 

It is of interest to monitor TW Hya with the spec- 
troastrometric technique over its full rotational periods. 
TW Hya is seen almost face-on. In the case that the 
alignment of the accretion columns to the position an- 
gle of the major axis of the disk is a sheer coincidence, 
we should see the locus of the astrometry rotating on 
the sky with the stellar rotation. TW Hya is the first 
young star-disk system where a planetary mass compan- 
ion might have been found (Setiawan et al. 2008), al- 
though its presence is questioned by the recent radial ve- 
locity measurements in the infrared (Huelamo et al. 2008; 
Figueira et al. 2010). The rotation period of TW Hya, 
inferred by the photometry and the line emission vari- 
ability, differs in the literature, ranging from 1.3 to 4.4 
days (Siwak et al. 2011; Mekkaden 1998; Dupree et al. 
2007; Batalha et al. 2002). Ill-determined rotational pe- 
riod critically hinders clear isolation of the radial velocity 
signal of the companion from the stellar activity. In the 
case that Br7 emission traces the accretion columns that 
are directly connected to the star, monitoring observa- 
tion with 2D spectroastrometry will help to visualize the 
stellar rotation and to provide the unambiguous rota- 
tional period. The present observations were performed 
on the two separate nights, on 2006 January 19 (MJD 
53754.20) and 2006 January 8 (MJD 53763.27). The in- 
terval was unfortunately too close to the multiples of the 
periods reported to date, and leaves little clues as of the 
stellar rotation. 
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Fig. 6. — Models of two-dimensional astrometry for disk rotation (left), accretion flow (middle), and disk winds (right). The upper 
panels show the expected centroid motion of Br7 emission. The co-rotational radius R c = 6.3 R* (Johns-Krull & Valenti 2001), where the 
disk rotation, accretion flow, or disk winds set off, is shown by dotted line. The filled circles at the center represent the actual size of the 
star (R„ = 0.96 Rq; Yang, Johns-Krull, & Valenti 2005). The lower panels show the observed line profile in dots, a Lorentzian fit to the 
observation in black line, and the line profiles that come out of the models in red. The assumed geometries of the ionized gas are illustrated 
in Figure 5. 
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